An Asymptotic Solution of the oo — d Hubbard Model 
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We present an asymptotically exact solution of the oo — d Hubbard model at 
a special interaction strength Ut corresponding to the strong-coupling Fermi-liquid 
fixed point. This solution is intimately related to the Toulouse limit of the single- 
impurity Kondo model and the symmetric Anderson model in its strong-coupling 
limit. 



PACS numbers: 71.30.+h, 75.10. Jm, 71.10.+X, 71.28.+d 



Introduction: Ever since the original introduction of the oo—d Hubbard model by Metzner 
and Vollhardt |jl]], and the recognition that the large-c? limit model has a fc-independent self- 
energy, the hope for an exact solution has been high. After a series of important analytical 
work 0- [^, and the recent surge in numerical simulations of the model, some of the 
outstanding physical processes described by the model have been gradually made clear. 
|T0| - |]14| . Yet, any exact solution is still out of sight. The bottleneck is: although the k- 



independent self-energy T,{uj) has reduced the problem to a + 1 dimensional one with a self- 
consistency, it is still highly non-trivial to calculate the Green's function of the corresponding 
"impurity" problem, which is crucial for the self-consistency to be complete. 

Here, we present a solution which is exact in the asymptotic limit, i.e. r ^ oo. The 
solution is obtained only at half-filling with a special interaction U. 

We first sketch the lines of reasoning reaching this solution: (i) The oo — d Hubbard 
model is mapped to a single impurity Anderson model with a self-consistent condition [§. 
The condition requires loosely speaking, that the local Green's function to have the same 
form as its surrounding conduction electron bath characterized by a hybridization function, 
(ii) As the interaction U is increased, the dynamics of the impurity model is governed 
by the "Kondo" physics, where there is an asymptotic point, called the Toulouse limit 
(TL) |l5l, at which the problem is exactly soluble. It will be shown that the corresponding 
symmetric Anderson Model also possesses such a limit, where the asymptotic behavior of the 
local Green's function can be obtained, (iii) With this Green's function, the self-consistent 
procedure can be carried out. It is found that, at the TL, the asymptotic behavior of the 
impurity Green's function can be made consistent with that of the hybridization function, 
thus the oo — d Hubbard model is solved asymptotically. 

Toulouse limit of the symmetric impurity Anderson model: It is well known that the 
single- impurity Kondo model has such an asymptotic limit: TL, where the conduction elec- 



trons become free |jT5||, [|T6[. Similar limit exists for the Anderson model. We first remind 
the readers how this limit is derived for the Kondo model emphasizing on the concept rather 
than the details. 



First, the Kondo Hamiltonian is separated into two parts, H = Hq + Hj, where 
-^0 = Y.k,a^kCj^„Ck,a + -fSzY.k,k'{Ck^^Ck',^ - Cfc^|Cfe',j) and Hj = J±T,k,k'{S'^C,^^^Ck,i + 
S~Cl^Ck,l)- Because of the spin-flip term Hj, a given spin sector of the system, say, the 
up spin, constantly undergoes transition between the eigen states of Hi = J2k ^kCl^-^Ck,^ + 
TSfc.fc'Cl.T^'fe'.T, and that of Hf = Y.k'^kCl^Ck,^ - xX^fc.fc'Cfc.tCfc'.T- The two eigen states 
have a relative phase shift, S = 2tan~^^^^^^, which is reflected in the propagator between 
the two states, {\e^^'^C^e~^^f^CJ\) ^ ^-(i-'^/'^)^. This is the origin of the non-Fermi hquid 
behavior of the X-ray edge problem |]l^, [|l^. When the phase shift takes a particular 
value, such that the exponent in the asymptotic Green's function becomes 2(1 — -)^ = 1, or 
pJz = |ct(7(-^), the Fermi- liquid (FL) behavior ^ is restored [jl9[, and the system appears 



to be free. This is the origin of the TL. The restoring of the FL behavior for an arbitrary 
phase shift 5 is a much more subtle result flrst pointed out by Anderson et al. [|1^ . 
The symmetric Anderson model 

H = Y. ^kClA^ + V Y.iClj. + h.c.) + Un\n{ - |(n( + nf ), (1) 

k^(7 k,(T 

does not have a spin-flip term explicitly. But via a discrete Hubbard- Stratonovich transfor- 



mation [20 1, we can write the partition function as. 



Wi=±i} 
where iE^e^'"^"^""^^ = e[^^"T"i+f(«T+"i)f^]^^ and coshX = e^^^^l Just like in the Kondo 
model, for a given spin sector, say the up spin, the system is alternating between eigen 
states of Hi = Ho + -tf\f^ (a = 1) and that of Hf = Ho - ^/|/t {(^ = -!)• Here 
Ho = Z]fc^fcC'fc,|C'fc,t + ^I]fe(C'fc,|/t + /|C'fc,t) • The phase shift of conduction electrons between 
these two eigen states is 5 = 2tan^^( J ^^ ) [ pT| . 

At 

For a given sequence of {cxi} = + + + H 1- + ++, where the flips are taking 

place longer than the relaxation time ||T^, |T6|, |21|, there is a well-deflned phase shift 
between the flipped states and the unflipped ones. This is reflected in the evolution operator, 
u{t) = {\e^^'^e~^f'^\). The reason that it is an evolution operator instead of a Green's 



function is because the spin-flip terms are absent here. This absence is compensated by the 
presence of the /—electrons in the Hamiltonian, which will contribute to an overall phase 
shift of TT, SO that (n(r)nt(r')) ^ (r - t')-^^-^^' still holds. 

To arrive at the TL point, the sum a is first regrouped into sections that contain equal 
number of spin-fiips, n. Within each section, the positions of the spin-fiips are summed, 
J2aia2(T3. ..(TL ~ ^n I ^ I ^ 5 ^^^ the avcragc is written as an exponential, UlS 



{-ir+Hn\^-I^\ 



(3) 



The partition function thus takes the form, 

~ ^qJo At Jo At "Jo Ar^ 

Once the TL is set, i.e. the coefficient in front of the logarithmic function becomes unity, 
the partition function is just the same as that of the following Hamiltonian: 

Ht = Y.^kClCk + VTj2iCld + h.c), (4) 

k k 

where Vr depends on the chemical potential of the Coulomb gas which is ignored so far, 
but will be introduced in the following calculations. Since Ht is free, the impurity Green's 
function < d(t)S{0) > can be easily obtained. 

Notice that the same Coulomb gas formula @ was derived long time ago by Hamann 



| 21| in his path integral approach. His tunneling configurations are the domain walls of the 
Ising variables here. That the symmetric Anderson model also has a TL is not surprising, 
considering the fact that via Schrieffer- Wolff transformation one can map it into a Kondo 
model. The subtlety is, this transformation allows one to go to a weak-coupling J of the 
Kondo, whereas the TL is a strong-coupling limit. The fact that the low-energy physics in 
the weak-coupling limit is controlled by the strong-coupling fixed point, keeps the physics 
of TL alive. 

Local Green's function of the symmetric impurity Anderson model: Comparing (1) to 
(4), the effect of the ex's sum is to push the effective / level to the Fermi point, i.e. zero, and 



to renormalize the coupling parameters such that the low energy behavior of the /-electrons 
is replaced by that of spinless (i-electrons. In fact, as pointed out by Anderson et al. ||16|| , 
(n°'(t)n^(0)) = {Sz{t)Sz{()))- This leads to a natural intuitive identification of d with / [^ . 
One can express the low frequency part of the Green's function as: 

{f{t)f\Q)),^^^^{d[t)sm, (5) 

where A/D is the spectral weight of the low frequency part. The weight of the high frequency 
part can be well approximated once the system is in the strong coupling regime, and the 
bare band shape is known. The sum of the two is unity obeying the sum rule. 

Another way of insuring the validity of the above relationship, is to use the slave-boson 
or slave- fermion decomposition scheme p3|, to separate /-electrons as fa = aba- + ad^bLa- To 



constrain the /'s in a singly occupied state, we restrict the bosons to J^a b^b^ = 1. It is then 
straight forwardly shown that Sz = ^(^j^t ~ ^I^i) ~ ^j^T ~ 1/2, S*"!" = 6|6j,S'^ = 6j6|, the 
usual Schwinger boson representation of the spin operators. If we neglect the doublon part, 
then we have < f^{t)f]{0) >= {a{t)a^{0)b^{t)b\{0)) ^ {a{t)a^{0)){b^{t)b\{0)) oc {b^{t)b\{0)), 
where one particle per site constraint has been enforced. On the other hand, we can write 
< (i(t)(i"l'(0) >=< S^{t)S^{0) > using the spinless fermion representation of the local spin 
operator. We then replace the spin operator by Schwinger bosons, such that, {d{t)d''{0)) =< 
b\{t)b^{t)b\{0)bi{0) >^< b^it)b\{0) >. Thus, we obtain eq.®. 

An exact asymptotic solution of the oo — d Hubbard model: As pointed out by Georges 
and Kotliar ^j, the most fruitful way to make use of a site-independent self-energy is to map 
the model to an Anderson impurity model plus a self-consistency. The impurity Lagrangian, 

£ = -jjdrdr'ClGo' C, + j drU{n^ - ^-){n^ - ^), (6) 

where Go is the Green's function with site 0, "the impurity site" removed. One first calculates 
G from £, and relates G to Go by, G~^{uj) = Gq'^^cu) — S(co'), thereby closing up the self- 
consistency. On a Bethe lattice, this relation is simplified to Go^^ = z — t'^G, The asymptotic 
form of Go can be obtained from its spectral density representation, Go(r) = J'^po{e)e~'''^de. 
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At zero temperature, Go{t) = ^sM _)_ 0{^). For oo — d Hubbard model, po{0) is pinned 0, 
because S(0) = 0, so that G{0) = Gq{0). On a Bethe lattice, the value of the pinning is ;^, 
where D is the radius of the semi-circle density of states. 

If we represent Go in terms of " integrated" conduction electrons in the Anderson Model, 
as we are allowed to do in the case of Bethe lattice, Gn ^ = z — V"^S— !— , we can transform 
(^) to (1) with a conduction electron density of states to be self-consistently determined. In 
this way, simply borrowing the last section can give the TL for the large-d Hubbard model. 
But here, we would like to derive this asymptotic limit in the Lagrangian formulation. The 
final results are the same. 

The partition function of the effective impurity version of the oo — d Hubbard model can 
be written as 

Z= Y: det{G^' + ^cr)det{G,'-^a), (7) 

Wi=±i} 

where the same discrete Hubbard-Stratonovich transformation is used. To obtain the 
Coulomb gas form, we make the following expansions of the determinants: 

detO^detO.^ = e-^'-(Goi;-)^+T.o(G^)...^ (g) 

where O^ = I + Gq-^ct- Since Go — * :;;^^, it is enough to keep the second order part 
only to obtain the asymptotic limit. Taking the trace in the imaginary time, Tr(GoO")^ = 
J dri J dT2Go{Ti — T2)(r{T2)Go{T2 — ri)a(ri) and using the fact Go(r) = — Go(— r), we have 

where /z is the chemical potential governed by the short time behavior of Go to be determined 



below. This long range 1-d Ising model can be mapped to a Coulomb gas model [0, p4 
The interaction strength in front of the logarithmics is 2( ^^^^ )^. When this is set to unity, 
we obtain the TL of the oo — d Hubbard model: J-^ = f^, where we have approximated 



A ~ vf/Ar. The low frequency part of the Green's function can be obtained from eq.(^): 

G(,),^ = ^/^ , « ^ (10) 



where Vt = ^, and T = 7cV^p{0) = ^V^. Since G'(O) = ^, we have A = ■^e~'^^', from 
which A can be determined, because /x is a function of A. 

Although the asymptotic solution G{t) = :^^ does not depend on A, neither does Ut, 
it is essential that A being finite [^ for the TL to exist as can be seen from eg . ([I0|) . To 



determine A, we have to consider the short time behavior or high energy part of Gq which is 
lattice dependent and cannot be calculated accurately within this framework. Fortunately, 
at least for the Bethe lattice, we can obtain an approximate value of A, which is finite at 
Ut, and thus proving the existence of the TL. 

At zero temperature, and for r > 0, Go{r) = J^ po{e)e~'''^de where po = —^IfnGo. The 
TL point is a strong coupling point. Our numerical experience with the Bethe lattice in this 



region shows that po consists of a (5— like peak around zq = vAD in addition to the finite 
part -i^ at z = 0. The former contribute to the short time behavior of Gq{t), the latter to 
the long time behavior. This has been derived analytically in Ref. [^. For completeness, 
we rederive it briefly here. On the Bethe lattice, Gq = ^_].2q - The G can be decomposed into 
a sum of low and high energy parts, G = Glow + Ghigh- Here Glow is the one obtained from 
eq.(^) and G^igh can be written as / ^^_^ '^ where ph only counts for the high energy part so 
that ph{t -^ 0) -> 0. For a small A, in the region A < z < \^AD, Ghigh{z) ~ -z/ ^^^^^ 



f2 ) 



and Giowiz) ^ ^. Thus, we have Goiz) ^ ,2^,^cyt^^/D ^^ere C = ej^andt = D/2. 
The position of the 6-\ike peak is found at, zq = ^i^^)^ with a weight of 27i+c)- ^O' ^^^ 
T being finite, we have: Go(t) ~ 2(1+0 ^~"^°^' Using (|) and (^, we obtain the chemical 
potential from the short time part of Gq: p ~ in+cF^"^^"^^" 

So far, everything is dependent on the cut-off At. There are certain arbitrariness in 
choosing this cut-off, which will affect the exact values of Ut and A. This is deeply rooted 
in the TL of the impurity models. Fortunately, the main result of this paper, i. e. the form 
of the asymptotic solution and the existence of this solution is independent of the choice of 
the cut-off. We adapt Hamann's choice of the cut-off Ar = ^ |21|]. With this choice and 
an estimate of C ^ ^7(1)^ = (^)^ we have, Ut = 2.QD and A = 0.2D. The result is 



very close to the second order perturbation calculation [|TT|, although as we pointed out, the 
exact number comparison is not very meaningful due to the cut-off dependence. 

For a hypercubic lattice, the kinetic energy is more spread out because of the Gaussian 
density of states in contrast to the semi-circular one. As a result, for the same value of 
interaction U, A is expected to be larger. Therefore, the TL will also exist there although 
the calculation in determine A is much more complicated. 

Discussion: (i) We want to emphasize that the TL point obtained from the Coulomb gas 
analogy is valid only asymptotically. In the Bethe Ansatz solution of the Anderson impurity 
model, there is no special value of U at which the rapidity equation becomes free. The same 
can be said about the Kondo model: even though there is a TL-like limit in the Bethe ansatz 



solution, it is different from the one derived by Toulouse [T^ or Anderson et al. |[T^. The 
difference lies in the phase shift. The Toulouse line in the Bethe ansatz solution intersects 
with the Wilson fixed point p2|. The phase shift of the Wilson fixed point is vr, the well- 
known unitary limit. The phase shift of the TL point is (1 — -^)7r, and is independent of 
the J±. In this sense, the usual TL does not show up in the Bethe ansatz for the Kondo 
problem either. 

(ii) The approximation involved here is to assume certain path being important, i.e. the 
spin-flips are separated long enough so that the system has a chance to relax. In this way, 
the concept of phase shift can still be used. In the Kondo regime, the interval of the adjacent 
spin-flip is about the order of the inverse Kondo temperature, which is long compare to the 
relaxation time. Therefore, this approximation becomes exact in the asymptotic sense. 

(iii) Within the above mentioned limitations, we have derived an exact asymptotic point 
of the oo — d Hubbard model, which corresponds to the strong-coupling FL fixed point. The 
form of the solution itself does not provide any new information about the model. This is 
due to a special feature of the large-d limit: the density of states is pinned on the metallic 
side. As long as the FL assumption is correct, one can arrive at this solution [H. But, the 



FL assumption can and does break down when the interaction strength is large enough [ITU 



The assumption has to be self-consistent to be valid. As stated in the beginning, this is the 



most difficult part of the oo — d Hubbard model. The contribution of this work is to show 
how the self-consistency can be achieved exactly at one point. To generalize to arbitrary 
interaction strength, we have also tried and failed to use another asymptotic approaches for 
this problem, namely, the conformal field theory |2^ . 



(iv) A more fruitful way of using the TL is to go away from half-filling, where much 
less has been understood. All the derivations used here are applicable except for one: the 
distinction between the high energy part and the low energy part spectral weight becomes 
blurred, thus making it more difficult to determine the value of A. 

(v) Through out the derivations, we have emphasized on the phase shift. Once the phase 
shift takes a special value, the model appears to be free. In this sense. X-ray edge problem 
also has a TL p8|, where the Anderson catastrophe seemingly goes away. But unlike the 



X-ray edge problem where the phase shift is well-defined between given initial and final 
states, for the single- impurity models and the oo — d Hubbard model, the phase shift is only 
meaningful in an intermediate step of a special treatment: the path-integral approach. In 
this special language, the system is phase-shifting back and forth until it is equilibrated. 



As creatively conjectured by Anderson et al. |]TB[, no matter what the phase shift is, as the 



system equilibrates, it will only take one value, the TL. This one statement turns out to 
involve all the machinery of many-body physics to justify and to modify [^. Thus, the 



essence of solving the model at the TL directly is to reach the strong coupling point without 
encountering the formidable task. The present work is just one more example in exploiting 
this limit Q. 
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